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Abstract 

A model for SU (n) string breaking on the lattice is formulated using strong 
coupling ideas. Although necessarily rather crude, the model gives an explicit 
picture of string breaking in the presence of dynamical quarks as a mixing phe- 
nomenon that involves the string state and a two-meson state. A careful analysis, 
within the model, of the Wilson loop shows that the evolution of the mixing an- 
gle as a function of separation may obscure the expected crossover effect. If a 
sufficiently extensive mixing region exists then an appropriate combination of 
transition amplitudes can help in revealing the effect. 

The sensitivity of the mixing region to the values of the meson energy and the 
dynamical quark mass is explored and an assessment made of the dectectibility 
of string breaking in a practical lattice simulation. 
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1 Introduction 



Currently much effort in lattice gauge theory simulations is being devoted to dynamical 
quarks |l|, For a review see S Giisken [0. For computational reasons the 

quarks incorporated in the model are still relatively massive. This makes it difficult 
to detect dramatic differences in the computed hadron spectrum relative to results 
obtained in the quenched approximation. In order to highlight the effect of dynamical 
quarks it is important to investigate phenomena that cannot occur in their absence. 
Gauge string breaking is an example of this type of event. Here we show explicitly 
that indeed string breaking does occur on the lattice and shows up, as expected, as a 
mixing phenomenon. 

The physical situation envisaged in lattice calculations comprises a static quark and 
anti-quark separated by a spatial distance of R lattice units (we mostly set the lattice 
spacing to unity). The two static particles may either support a gauge string that runs 
between them or separately bind a quark and an anti-quark to create a two-meson 
state. When the string is stretched sufficiently, its energy coincides with the energy of 
the two meson state. In the neighbourhood of this critical separation the two different 
physical states mix permitting the occurence of transitions between them. 
There are two related issues of importance, namely the size of the range in R over which 
there is substantial mixing and the magnitude of the energy split induced by mixing. 
As fractional effects we find that they are proportional to one another and therefore 
of essentially the same size. We show in terms of our model that, for the lowest quark 
masses and for meson excitations that are not too great, the mixing effects of string 
breaking may well be detectable. If the meson masses are too great and/or the quarks 
too heavy then the phenomemnon will be undetectable on the lattice. 
Because it is based on strong coupling ideas the results of our model are provisional 
and require interpretation if they are to be applied to actual simulations. Nevertheless 
the model is very suggestive and identifies crucial aspects of the mixing phenomenon 
to which the results are sensitive. 

2 Strong Coupling Model 

We formulate the model by first recalling the well known rules for evaluating simple 
graphs in the strong coupling limit of SU{n) gauge theory ^. They are as follows: 
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1. A factor of 



for each plaquette. 
2. A factor of 



2k 



2n2 



1 + 7.e 



2 

for each Wilson quark hne in the direction of the unit vector e . Here k is the 
standard quark hopping parameter. 

3. A factor of (—1) for each internal quark loop. 

4. A factor of 1/n for each internal quark loop. 

5. A trace over the spin matrix factors for each internal quark loop. 

2.1 Simple String Model 

In the absence of dynamical quarks our model for the correlation function of the string 
of length R over an imaginary time interval T is the standard R x T Wilson loop. 
In leading strong coupling approximation the above rules give for this string-string 
propagator 



Gss{T) 



R 



T 



(1) 



2n^ J 

This tells us immediately that the energy of the string state is V{R) where 

V{R) = aR , (2) 
and a is the dimensionless string tension given by 

--log(^) . (3) 

Note that we have omitted the spin degrees of freedom of the static quarks. They play 
no role in the model. 
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2.2 Model for Mesons 



The simple model for mesons suggested by the strong coupling limit, is one in which a 
light quark propagates along a static (anti-)quark line. See Fig 1 . Using the rules, we 
find that this meson propagator has the structure 

9{T) = ^^{2k)' . (4) 

Of course we retain the spin degrees of freedom of the light quarks. They do play a 
role in the model. 

The propagator for two mesons moving independently each bound to its static quark 
is 

= (ii^)"' <8 {^f ((2K)f . (5) 

In fact, as will become clear below, only a particular combination of quark-anti-quark 
spin wave functions enters the calculation. It is obtained by completing the light quark 
loop with the matrices 



l±7i 
2 ' 

and including a factor of (—1) . This two- meson propagator is 



(6) 



'MM 



(T) = -Tr (i^) g«"(r) (l±2l) ,(2)p.) . ]. ((2«)f . (7) 



Finally we will give the light quarks, when bound in a meson, a hopping parameter 
k' that may be distinct from the value k that holds elsewhere in the diagrams of the 
model. Eqd^) should be modified therefore by the replacement k, —>■ k' . The essential 
point is that we are treating the energy, Em, of the mesons as a parameter. This is 
reasonable since the energy of the mesons is an important number in the string breaking 
process and a simple version of the hopping parameter expansion cannot refiect the 
full dynamics of the meson bound state. We asume then 

(2/«')' = e-^" , (8) 

where Em is the combined energy of the two bound mesons. 



3 String Breaking 

Our model for string breaking involves summing over planar graphs that incorporate 
transitions between the string and meson states described above. We arrive at the 
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rules for computing these graphs by considering the case of string-to-string propagation 
represented by Fig 2. Prom the rules we see that the corresponding contribution to the 
Green's function is 



2n2 



./\2 



T2 



where 



t = Tr 



{2k') 

\2Kf 

1 + 7o\ /I + 71 



2v? 



R 



R' 

I 

t 
11 



/\2 



{2k') 



{2k) 



T4 



2n2 



7o 



71 



1 

'2 



(9) 



(10) 



The factors involving k (as opposed to k') are associated with the vertical sides of 
the quark loops. Note that the trace of the internal quark loop is equivalent to the 
quark-anti-quark spin projection mentioned previously. 
We reinterpret the structure of such diagrams as follows: 



1. The factor 



2v? 



R 



propagates the string by one time step. 
2. The factor 

{2k 



/\2 



propagates the two-meson state by one time step. 
3. The factor 



1 



{2k) 



R 



is associated with the transition from string to two-meson state and vice-versa. 

At any stage we may view the system as being in either a string or a two-meson state. 
To describe the transition from an initial time to time T we need a 2 x 2 matrix of 
transition amplitudes 



G{T) 



Gss{T) Gsm{T) 
Gms { T) Gm m{T) 



If for convenience, we introduce the parameters 

R 



a 



2-0? 



b = {2k') 



{2k) 



R 



(11) 



(12) 
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then the above stepping procedure can be represented by 

^^^^^^-^\Gms{T) Gmm{T) ) ' 
where the matrix A is given by 

A - ( 

^-[bc b ) ■ 

If for definiteness we set 

G(0) = 1 , 

then 

G{T)^{Af . 

The matrix A is not symmetrical but can be expressed in terms of a 
B in the form 

A^DBD-^ , 

where D is the diagonal martix 

D-( \ 

[ Vb) ' 

and 

^ _ ( c \ 

\ \/ab c b J 

The eigenvalues of B are 

i|(a + 6)±^(a-6)2 + 4a6c2| . 
The corresponding eigenvectors are 
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_ / cos 6 \ _ I —sin 6 

^+ ~ I sin^ I ' ^" ~ I cos^ 



-(a -b) + ^{a - bf + 4a6c2 
tan 6 — 



where 



2\/ ab c 

If we construct the othogonal matrix 

Q _ ( cos 9 — sin ^ 
I sin^ cos^ 



then 

B = OAQ-^ (24) 

where 



It follows immediately that for any T 



A^^DOl^ ^+ j O-'D-' . (26) 

This shows clearly that the propagating eigenstates are mixtures of the string and 
two-meson states determined by the mixing angle 9 . The corresponding eigen-energies 
are 

£;± = -logA± . (27) 
Note that is the lower of the two energies. 

4 Wilson Loop 

We can now use the above results to elucidate the behaviour of the Wilson loop provided 
we pay attention to end conditions. The Wilson loop begins with a factor a . The sum 
of all such diagrams of length T are obtained as the components of 



iT-l ( « 





(28) 



The Wilson loop itself is obtained as 



W{R, T) = GssiT) = (1, 0)DOA^-'O-'D-' j . (29) 

That is 

W{R,T) = a {cos'' e Xl-^ + sin'' e X'r^) . (30) 
The Wilson loop, therefore, sees both exponentials 

e-^±^ = A^ . (31) 

In principle, the the lower energy exponential should dominate the asymptotic be- 
haviour. In practice, what is observed will be influenced by the behaviour of the 
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coefficients of the exponential contributions through their dependence on the mixing 
angle. 

Noting that 

-{a - h) + J (a - by- + l«&c2 

8111^9 = — ^ , (32) 

2^(a - 6)2 + 4a6c2 

we see that when V{R) « Em then a » b and we have 

^~0, A+~a, A_~6, (33) 

with the result that Ej^ ~ ^(-R) and ~ Em ■ Under these circumstances the 
coupling to the state with energy E^ will vanish and only the exponential associated 
with £■+ will be observed as expected. When R is chosen so that V{R) >> Em then 
b » a and we have 

, X+^b , A_~a , (34) 

with the result that E^ ~ Em and E_ ~ V{R) . Under these circumstances the 
coupling to the state with energy E^ will vanish and only the exponential associated 
with will be observed. That is, both above and below the crossover point, only the 
original string behaviour 

g-y{R)r 

will be observed. The movement of the mixing angle therefore obscures the crossover 
phenomenon in which the string energy is expected to be bounded by the energy of 
the two-meson state. As is clear from the above analysis this does not happen. 
To see mixing effects directly in the Wilson loop we must fix R near the critical crossover 
value Rc, for which V{R) = Em and hence a — b . That is 

Rc-^ . (35) 

At this point 9 — tt/A and 

WiR,T) = "-{xl-' + Xl-') , (36) 

so the two exponentials appear with equal strength. In order to separate them in a 
measurement it must be possible to measure their difference. We find 

A± = 6(l±c) . (37) 
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So approximately for small c 



E± = -logAi ~Em(1tc) , (38) 

that is 

Em cfRc ^ ^ 

We can define the mixing region as the range of R for which both sin^ and cos^ 6 are 
significantly different from unity. The range of the mixing region in R can be estimated 
as 

(40) 



Rc \2R d9 
It is closely related to the mixing energy we find 

AR AE 

IX- 



R—Rc 



(41) 



Rc Em 

We estimate the numerical significance of these results below. 

5 String-Meson Transition 

The string-meson transition amplitude can computed from eq(p8|) as 

Gms{T) = {QA)a(^1^ . (42) 

This gives 

GMsiT) = Vah sine cos 6 (A^^^ - A^"^) . (43) 

Because of the presence of the factor sin 6 cos 6 this transition amplitude vanishes out- 
side the mixing region. In this region, however, it is possible to detect the presence of 
the dominant energy provided the two energies are sufficiently well separated. However 
if the two eigen-energies are only weakly separated then the amplitude will experience 
a further suppression because the two contributions have opposite signs that give rise 
to a cancellation. This is consistent with the results of the previous section on the 
Wilson loop and fits in with the obvious idea that if the primitive transition process 
vanishes then the transition amplitude itself vanishes. 
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6 Meson-Meson Transition 



The meson-meson transition involves a sum over amplitudes that begin and end with 
a factor b . They can be evaluated as 

gMM{T) = {0,l)A^l^ . (44) 

The result is 

gMj^{T) = b (sin^ 9 A^-i + cos^ 9 \^_-') (45) 

Clearly this amplitude shows properties similar to that for string-string transitions 
in that on both sides of the mixing region in R only the meson energy is detected 
in the T-dependence of the correlator. However the result suggests that when R is 
within the mixing region, the combination b QssiT) + a Qmm{T) will show a reduced 
influence from the mixing angle and will allow the best chance of tracking both and 
individually from the resulting T-dependence. Of course this requires favourable 
circumstances in which the mixing region is sufficiently large and the energy separation 
of the mixed states is sufficiently great. 

7 Interpretation of Results 

In order to use the results of our analysis wc will make the assumption that the pa- 
rameters of the model can be interpreted directly in terms of the appropriate physical 
parameters of an actual simulation rather than as bare lattice parameters. We have 
already anticipated this approach when we introduced k' as a way of parametrising the 
meson energy Em ■ In the same spirit wc will interpret the other hopping parameter 
n in terms of the light quark mass, m^, and require 

2k = 6-"^^ . (46) 

We therefore have the identifications 

a = e-"^ , b = e-^« , c = ^e""*^^ . (47) 
With these identifications we then see that the fractional energy shift due to mixing is 

Em V 3 Em 



10 



niq (GeV) 


AE/Em 




.1 


.7077 


18.3 


.2 


.4089 


10.6 


.3 


.2363 


6.1 


.4 


.1365 


3.5 


.5 


.0789 


2.0 



Table 1: These results are for a"^ = 1.5 GeV, a^/^ = 0.427 GeV, Em = 1-0 GeV. The 
critical distance is Rc — 8.22 lattice spacings. 



(GeV) 


AE/Em 


AR 


.1 


.2044 


10.6 


.2 


.0682 


3.5 


.3 


.02295 


1.2 


.4 


.0076 


0.4 


.5 


.0025 


0.1 



Table 2: These results are for a"! = 1.5 GeV, a^/^ = 0.427 GeV, Em = 2.0 GeV. The 
critical distance is Rc — 16.45 lattice spacings. 

and the mixing range in R is 




AR = n\ - . (50) 



We have now specialised to SU (3) . 

If we recall that Rc — Em /(J then we see from this formula, through its dependence on 
Rci that the mixing range is very sensitive to both the meson energy, Em and to the 
quark mass, ruq . In fact 

(7 

In order to obtain numerical estimates we consider a range of reasonable values appro- 
priate to practical simulations for these parameters. The results are shown in Tables 
1& 2 . These results show directly the sensitivity of the physical picture to the energy 
of the two-meson state and to the quark mass. For an inverse lattice spacing of 1.5 
GeV and a two-meson energy of 1 GeV, we obtain a critical separation Rc of roughly 8 
lattice spacings. At the rather light value of 100 MeV for the quark mass we see that 
the mixing region is large and covers a range of more than twice this critical distance. 
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For such circumstances we might expect to see the effects of mixing clearly. As the 
quark mass increases in size the mixing region shrinks rapidly so that by about 400 
MeV it has reduced to a barely observable range of just over 1 lattice spacing either 
side of the critical separation. When the two-meson energy is increased to 2 GeV only 
the lightest quark mass of 100 MeV yields a mixing region of any size. At the lower 
two-meson energy the the results suggest that there would be a reasonably detectable 
energy difference for the lower quark masses. At the larger two-meson energy only the 
lightest quark mass yields a substantial mass difference in the mixing region. 

8 Conclusions 

We have formulated a simple model of string breaking in terms of strong coupling 
ideas. The model shows explicitly string breaking as a mixing phenomenon between 
string and two-meson states. Because of the nature of the dependence of the Wilson 
loop on the angle of mixing between these two states, we do not expect to see the 
crossover phenomenon for the string tension in this amplitude. Instead the original 
string tension should be observed on both sides of the mixing region in i? . In a similar 
way the meson-meson correlator will show only the meson energy on both sides of the 
the mixing region. The string-meson transition amplitude vanishes outside the mixing 
region and suffers a suppression inside the region because the eigen-energy exponentials 
enter with opposite signs. This amphtude should therefore show at best only a weak 
signal in a simulation unless the mixing region is sufficiently large and the eigen-energies 
sufficiently distinct. The model suggests that a particular combination of the string- 
string and meson-meson amplitudes would show a reduced dependence on the mixing 
angle and give the best chance of tracking the two eigen-energy exponentials through 
the mixing region. In a practical simulation this combination would have to be sought 
empirically. 

To observe string breaking it is necessary to be able to measure amplitudes in the 
mixing region with sufficient accuracy. For this to be possible a simple interpretation 
of the model in terms of realistic simulation parameters suggests that it will only be 
possible if the meson energy is not too high and the quark mass is sufficiently low. More 
precisely we would expect that at the critical value of the static quark separation, i?c, 
we have V{Rc) < 1.0 GeV and niq < 200 MeV. The sensitivity of the mixing region to 
the value of the light quark mass, predicted by the model, may yield interesting results 
on string breaking if exploited systematically in practical simulations. 
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Figure 1: A dynamical quark (light line) bound to a static quark (heavy line). There 
is closure at the ends of the lines and an implied trace for the gauge degrees of freedom 
but not for the light quark spin degrees of freedom. 
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Figure 2: Wilson loop (heavy line) containing internal quark loops (light lines). 
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